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Abstract. As it is well known for stiff equations of state an existence of neutron stars with masses ^ 2 Mq is 
possible. Especially interesting possibility is opened if the equation of state based on the Skyrme theory is realized 
in nature. This equation of state was proposed recently by Ouyed and Butler. We discuss different channels of 
formation of massive rapidly rotating neutron stars. We use a population synthesis code to estimate numbers 
of massive neutron stars on different evolutionary stages. A neutron star increases its mass by accretion from 
a secondary companion. Significant growth of a neutron star mass due to accretion is possible only for certain 
values of initial parameters of the binary. In this paper we show that significant part of massive neutron stars 
with M ^ 2 Mq can be observed as millisecond radio pulsars, as X-ray sources in pair with white dwarfs, and as 
accreting neutron stars with very low magnetic fields. 
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1. Introduction 

Mass is one of the key parameter for neutron star (NS) 
physics and astrophysics. It can be measured with high 
precision in binary radio pulsar systems. Up to very re- 
cent time obtained results fell in a very narrow region 
1.35-1.45 Mq Ixhorsett & ChakrabjJty 1999). These val- 
ues lie very close to the Chandrasekhar limit for white 
dwarfs. Thus, for years M = 1.4 Mq was considered to be 
a standard value of a NS mass. Recently the range widened 
towards lower mass es thanks to the disc overy of the double 
pulsar J0737-3039 l|Burgav et a l."2003'). One of the NSs in 
this system has M = 1.25 Mq iFyjie et al. 2004). Also the 
mass range expanded towards higher masses, though this 
result is less certain. There is only one NS in a binary ra- 
dio pulsar system with mass significantly higher than the 
canonical value 1.4 A/q. It is the pulsar J 0751-I-1807 with 
the m ass 2.llg g (95% confidence level) ijNice &: SplaveJ 
iool. All others are consistent with the standard mass 
value on the 1-2-a level. However, small number of mas- 
sive radio pulsars can be a result of a selection effect (s). 

There are reasons to suspect an existence of signif- 
icant number of NSs with higher masses. Evidence for 
such a proposal comes both from theory and observations. 
Calculations of cooling curves of NSs suggest that some of 
these objects might be more massive tha n known sources 
in ra dio pulsar systems (see for example, * Kaminker et alJ 
l200l|) with M up to 1.8 AIq and probably more. Modeling 



of supernova (SN) explosions also suggest the existence of 
NSs with higher masses (^ooslcv ct al. 2002). Still mod- 
els of NS thermal history and SN explosions do not requier 
masses M > 2 Mq, but there are observational indications 
for their existence. 

Obscrvationally high masses of NSs are mainly sup- 
ported by data on X-ray binaries (we do not discuss 
here data based on quasi-periodic oscillations measure- 
ments as they are model dependent). Estimates for sev- 
eral system s give very high valu es: 1.8-2.4 Mq for 
Vela X-li llQuaintrell et all l2003l) . 2.4± 0.27 M,7^ for 



4U 1700-37 Iciark et all l2002t see also jHeineke et al 



l2003t Ivan KerkwiikI 1200^ . Very recentlv Ishahbaz et al 



(|2004|) presented observations of a low-mass X-ray binary 
2S 0921-630/V395 Car for which 1-a mass range for the 
compact object is 2-4.3 Mq. Still it is necessary to note 
that such measurement are not as precise as the radio pul- 
sar ones (for example, at the 3-fT level the mass of the NS 
in Vela X-1 is still compatible with the standard value). 

The existence of NSs with M ~ 2-2.4 Mq is not in 
contradiction with the present day theory of NS interiors. 
There are several models with stiff equation of state (EOS) 
which allow an existence of NSs with masses ^ 2 Mq (see 
a review and references in Hacnsel 200^. Here we will 
focus on so called skyrmion stars (SkyS) as they are ex- 
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This range is ba sed on the two estimates given 
iQuaintreU et al.ll2003^ : 1.88 ±0.13 and 2.27±0.17Mq. 



2 



S.B. Popov and M.E. Prokhorov: Formation of massive skyrmion stars 



pected to be a kind of NSs with the highest value of max- 
imum mass (Mmax)- 

In 1999 Ouye d and Butler d iscussed an EOS based 
on the model of ISkvrmel l|l962l) . A NS with such EOS 
has Mniax=2.95 Mq even for a non-rotating configura- 
tion. Usualy maximum rotation can increase the limit 
by ^ 15-20%. Rapidly rotating SkyS were discussed by 
Ouyed (2002, 2004), and for this case Mmax=3.45 Mq 
and i? = 23 km (this model also has relatively large radii 
of NSs). Such model is very interesting from the astrophys- 
ical point of view, and it is important to discuss scenarios 
of formation of compact objects with such high masses. 
Our goal in this note is to pick out evolutionary tracks 
of binary systems which can lead to the formation of NSs 
with high masses, and to discuss possible observational 
appearences of such systems and their relative and abso- 
lute numbers in the Galaxy. As we do not use explicitly 
any EOS in our calculations, then our results can be ap- 
plied to other stiff equation of state and even to low-mass 
black holes (BHs). 

In the next section we discuss evolutionary paths at the 
end of which a massive NS can be formed. Then we give 
an estimate of the number of massive NSs in the Galaxy. 
Finally we discuss our results and propose systems which 
are more favorable to host massive NSs. 



2. Possible channels of massive neutron star 
formation 

As mass determination for NSs is possible only in binary 
systems^ we focus on potentially observable stages of evo- 
lution of binary systems in which a massive NS can form. 
Below we discuss possible ways of massive NS formation. 

Since we are mostly interested in compact objects with 
rapid rotation (because they can have higher maximum 
masses) it is necessary to follow evolution in a binary as 
such objects cannot form from single stars teeger et alJ 
l2003j) . so its necessary to study evolution of close binary 
systems. Except evolutionary tracks which lead to a for- 
mation of a massive NS in a binary we follow paths at 
the end of which an isolated massive NS can form. An 
appearence of a rapidly rotating single massive NS due 
to a binary evolution can be a result of a coalescence of 
two compact objects (NSs or white dwarfs — WDs), or a 
result of a more slowly merging process in which a normal 
star is involved, or a result of an evaporation of a low-mass 
secondary companion by an active pulsar. At some stages 
during its evolution a binary which finally is going to pro- 
duce an object of our interest can be observed as an X-ray 
source, that is why it is important to select evolutionary 
paths also for them. 

The main output of a collapse of cores of massive 
stars are NSs with M ~ 1.2-1.5 Affy. This conclusi on is 
supported both observationally I van Kerkwiikll20f)3) and 



^ Note, that in principle there is a possibility to determine 
an isolated NS mass by microlensing effects, however, we do 
not touch this issue here. 



theoretically llTimmes et al.ll996HFrver fc Kalogeral200ll: 
IWooslev et al .I2OO2I) . Numerical models of collapse are not 
as precise as necessary to determine the exact shape of a 
NS mass spectrum (for example the amount of fallback 
is not well known), however, calculations show that the 
formation of NSs with high masses is not favourable and 
most of them should have M ^ 1.3-1.4 Mq. 

A discovery of a NS with M 1.8 A^q should mean 
that the mass was increased after formation of the com- 
pact object during its evolution (if the mass is significantly 
higher than 1.8 Mq then such a conclusion seems to be 
inevitable). Based on this proposition we call below as 
massive NSs with M > 1.8 Mq. 

A NS can increase its mass due to fallback, coalescence 
with another NS, or accretion from a secondary compan- 
ion. As we note above, the first way is not well studied, 
and we do not discuss it be low. Oppositely coal escence of 
NSs is well understood fsee lRosswog et al.l2003l and refer- 
ences therein). The rate of NSs coalescence in the Galaxy 
is about 1 per lO"' yrs. As a result a rapidly rotating mas- 
sive isolated NS (or a BH) can form. This way of evolution 
also will not be discussed below. In the following only bi- 
nary evolution of a NS in pair with a normal star or a WD 
will be studied. 

At first for an illustration let us assume an isotropic 
collapse, ie. zero kick. Such an assumption is not realis- 
tic as most part of NS - nearly all radio pulsars - ob- 
tain at birth high addition al velocity ^100-1000 km s~^ 
l|Arzoumanian et alJl200^ . However it is much easier to 
understand main processes in a binary evolution if one 
neglects kick. In addition, if a binary was not unbounded 
after a SN explosion then an orbital eccentricity quickly 
decays after the secondary fills its Roche lobe. So, if at the 
moment we are not interested in the question of the bi- 
nary survival then it is possible to neglect kick to simplify 
the explanation. 

Let us start with a qualitative discussion (below in 
sec. 2.1 a more detailed consideration is given). The 
most obvious channel to form a rapidly rotating mas- 
sive NS is an evolution in a low-mass or intermediate 
mass binary (see, for ex ample, recent calculations by 
'Podsiadlo wski et a l."2002^. This path includes, for exam- 
ple, millisecond pulsars (however it is not the only possible 
output). 

As we are interested here in systems with high mass 
ratio (a massive primary produces a NS and the secondary 
star has a low mass) it is necessary to consider three differ- 
ent situations after the NS formation when the secondary 
fills its Roche lobe: i.) a normal star can fill its Roche lobe 
without a common envelope formation; ii.) a normal star 
can fill its Roche lobe with a common envelope formation; 
Hi.) a WD fills its Roche lobe. 

To fill the Roche lobe a normal secondary star has to 
evolve further the main sequence stage. During its evolu- 
tion prior to the Roche lobe overfiow the mass of the star 
is nearly constant (see detailed tracks below). A common 
envelope is not formed if the normal star is not signifi- 
cantly heavier than the NS. In this regime mass is not lost 
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from the binary system. For more massive secondaries for- 
mation of a common envelope is inevitable, mass transfer 
is unstable. In this regime significant fraction of the mass 
flow is lost from the system, so the mass of the NS grows 
less effectively. It is only partly compensated by higher 
mass of the donor. 

After the common envelope stage an orbital separation 
becomes smaller, so later on even a degenerated core of the 
secondary - a WD - can fill the Roche lobe. 

2.1. Evolutionary tracks 

For our calculations we use the "Scenario Machine" code 
developed at the Sternberg Astronomical Institute. ^ 
De scription of most of pa rameters of the code can be found 
in ijLipunov et alJll99 6^. Below we mention those which 
are the most important for us here: 

— All NSs are born with M = 1.4 Mq. 

— At the common envelope stage a hypercritical accre- 
tion (with M much larger than the Eddington value) 
is possible. 

— During accretion the magnetic field of a NS decays 
down to the value which cannot prevent rapid (maxi- 
mum) rotation of the NS. 

— Oppenheimer-Volkoff mass of a rapidly rotating NS 
(the critical mass of a BH formation ) is assumed to be 
3.45 Mq according to lOuvedl l)2004j) . 

For zero kicks we distinguish two groups of tracks 
which produce massive NSs. A typical track from the first 
group has initial value of the semimajor axis a = 290 Rq 
and star masses Afi = 10.5 Mq, Ma = 2 Mq (fig. 1 left)^. 
After the massive component leaves the main sequence it 
expands and fills its Roche lobe. As a result the common 
envelope stage sets on. During this stage the orbit shrinks 
by more than an order of magnitude, and the primary 
looses about 3/4 of its mass and becomes a low-mass he- 
lium SN progenitor. After the SN explosion the orbit has 
low eccentricity and a ~ 7-8 Rq. Mass of the secondary 
is not changed during these stages of the evolution. 

Till the secondary fills its Roche lobe the NS is at the 
stage s of ejector and propeller (see for example iLiounovl 
ll992l for stages descriptions). During these stages the mag- 
netic field is assumed to be const ant. Stage durations can 
be found in ((Lipunov et al.lll99(t ^. 

After the secondary fills the Roche lobe the NS starts 
to accrete. At that moment the mass ratio is about 0.7 
(the NS is lighter) and a mass transfer is stable with nearly 

^ Online materials are available at 

fhttp:/ /xray.sai. msu.ru/sciwork/scenario. html' and 
[http:/ /x ray.sai.msu.ru~~mystery/articles/review/ 

Colored version of the figure in 
high resolution is avalable on the Web: 
|http:/ / xray.sai.msu.ru / ~polar/htmI / publi cations /ouyed/| 

The subsonic propeller stage is not taken into account as 
for binaries with big accretion rates it is very short. 



zero mass loss from the system. Up to equalizing of com- 
ponents masses matter transfer goes on a thermal time 
scale, after equalizing - on a nuclear. Process of accre- 
tion can be stopped because of a switching on of a mil- 
lisecond radio pulsar. It happens when the donor's mass 
is '~ 0.1 Mq. The remnant of the secondary companion 
then can be evaporated completely, while the evaporation 
is proceeding the systems looks like the famous "Black 
widow" pulsar 1957-1-20 (and its twin PSR J2051-0827). 
If accretion is not stopped then it continues till a planet- 
like (Jupiter mass) companion remains. As we see the final 
stage of such an evolution is a "single" massive rapidly ro- 
tating NS. In both cases the final mass of a NS can reach 
3.2-3.3 Mq. We can observe such a system at the stage 
of accretion which lasts 90% of the evolution. Masses of 
NSs in these accreting systems can be in the range from 
the initial mass (1.4 Mq in our case) up to 3.2-3.3 Mq. 
Orbits can be relatively wide. 

The described evolutionary channel appears to be nar- 
row in a sense that small changes in the initial conditions 
do not allow a massive NS formation. Also uncertain pa- 
rameters of the common envelope stage can significantly 
influence this path. 

Ranges of initial parameters of evolutionary tracks 
from the second group are given in the table 1. We give 
maximal and minimal values for two types of tracks (2a 
and 2b) which differ by the final stages of evolution. The 
orbital period, Porb, is given in the table 1 just for an illus- 
tration. In our calculations we use masses and semimajor 
axes. So the values of Porb given in the table are simply 
calculated using maximum masses and minimum semima- 
jor axes for the shortest period, and, oppositely, minimum 
masses and maximum axes for the longest period. Due to 
that ranges for Porb for tracks 2a and 2b intersect. 

A typical representative of the 2a subgroup has the 
following initial parameters: a — 300 Pq, Mi — 12 Mq, 
M2 — 4Mq. The main difference form the first group of 
tracks is a more massive secondary companion. Because 
of that the common envelope during the first mass trans- 
fer is less effective, and after a SN a system with a = 
170 Pq and low eccentricity is formed (the mass of the 
secondary is not changed). Later the secondary fills the 
Roche lobe. Mass ratio is high, mass transfer is unstable 
and the common envelope forms. At the end of the com- 
mon envelope stage the secondary becomes a WD with 
M ~ 0.8 Mq, and the orbital separation diminishes down 
to 5 Pq . During the common envelope stage the NS in- 
creases its mass up to ^ 2.3 Mq (for more massive donors 
mass loss from the system is more effective, so in such 
cases the NS mass can be lower: ^ 1.9 Mq). 

After the formation of a binary consisting of a NS and 
a WD the evolution in the second group can take one of 
two different paths. For some tracks (2a) from the sec- 
ond group the time of rapprochement of the components 
due to gravitational wave emission is too long, so there 
is no Roche lobe overflow. Systems with smaller orbital 
separation have enough time to approach to each other 
close enough for the beginning of WD overflow. This sit- 
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Table 1. Parameters for tracks from the second group 



parameter 


min 


max 


width 


Track 2a 
a 

Ml 
Mi 

-Porb 


279% 
IO.SMq 
3.9Mq 
123"^ 


670i?Q 
12.8M0 
6.7M0 
537"^ 


0.20 

0.054 

0.13 


Track 2b 
a 

Ml 

Ma 

P K (*' 

orb 


135% 
1O.3M0 
3.9M0 

4ld 


279% 
12.4M0 
7.4M0 
144d 


0.17 

0.046 

0.15 


-Porb is i 

the text. 


jiven just as an illustration, see 



nation corresponds to the initial parameters a = 200 Rq , 
Ml = 12M0, M2 = 4M0 (track 2b in the table 1). 

The main difference between tracks 2a and 2b is 
smaller orbital separation in the latter case. Track 2b is 
similar to the one on the right panel of fig. 1 , but after the 
common envelope semiaxis of the system is just ~ 3-Rq. 
A WD has enough time to fill the Roche lobe and com- 
pletely transfer its mass to the NS. At the end we have a 
single rapidly rotating NS. The NS mass for this case is 
increased up to ^ 3 Mq. Stages with a WD are shown in 
the box as they distinguish the track 2b from 2a. 

When a semimajor axis is larger than a ~ 670 Rq 
the second common envelope results in NS-star merging, 
so the Thorne-Zytkow object is formed. Its evolutionary 
path is not very clear. A formation of a massive NS and 
a formation of a BH are both possible. We do not include 
this possibility into our calculation. 

2.2. Evolutionary tracks with kicks 

Above we discuss two families of tracks with zero kicks 
which result in massive NSs formation. However, it is nec- 
essary to include kicks as they are a general property of a 
NS formation. A kick can change orbital parameters after 
a SN explosion, it can even make the system unbounded. 
If after a SN (and after a brief period of circularization 
of an orbit) we obtain in our calculations a system with 
parameters in the range which was obtained above for the 
zero kick, then the following history of the system should 
be the same as described in sec. 2.1. 

An additional velocity which a NS obtains at birth can 
change the range of initial parameters that are necessary 
for a massive NS formation. Especially it is important to 
estimate if ranges for Mi, A'l2 and a are changed signifi- 
cantly or not. As a kick velocity and a NS mass in our cal- 
culations are assumed to be independent on a mass of an 
exploding star (see below sec. 4.1) we do not expect that 
a range of masses of primaries should be modified. The 
same can be said about the range of initial masses of sec- 
ondaries because these stars do not suffer any important 
evolutionary changes before a SN expolsion. Since a kick 




Fig. 1. Evolutionary tracks for massive NS formation. 
In the left panel we show a typical track from the first 
group. The first mass transfer (from the primary) results 
in a common envelope formation due to high mass ra- 
tio. Accretion onto a NS from the secondary companion 
proceeds stably without a common envelope. In the right 
panel we show an evolutionary path of a system from the 
second group. This track differs by a higher mass of a 
secondary companion. Because of this difference the first 
mass transfer goes on without a common envelope. A NS 
gathers an additional mass during one or two episodes 
of accretion. If the orbital separation is not very large 
(~ 200 i?Q, see text) then at first the NS accretes from a 
normal secondary filling its Roche lobe, and then from a 
WD (this stage is shown in the dashed frame). For wider 
systems the evolution stops after the mass transfer from 
the normal secondary star (ie. before the frame). On the 
left from each track we indicate evolutionary stage (in no- 
tation from Lipunov, Postnov & Prokhorov 1996) and NS 
masses. 

can dramatically change the orbital parameters the situ- 
ation is different for the initial orbital separation range. 
For example, with a kick systems wider then the ones dis- 
cussed in sec. 2.1 can still form massive NSs. 

In the next section we present results of our calcu- 
lations of population synthesis of massive NSs for both 
scenarios. 
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Table 2. Fractions of massive NSs at different stages 



Stage with kick without kick 

Ejector 032 0.39 

Propeller + Georotator 0.02 0.08 

Accretor 0.66 0.53 



Hypercritical stages 5 ■ 10 



3. Estimate of observable number of massive 
neutron stars in the Galaxy 

To estimate the number of massive NSs in the Milky Way 
we run several sets of population synthesis calculations for 
the ranges of initial parameters which correspond to the 
two groups of tracks described above. Each run includes 
calculations of 10^ individual binary evolutionary tracks. 

We run the model for zero kick velocities and for non- 
zero ones. For the latter ca se we use the distribution simi- 
lar to the one suggested in ((Arzoumanian et al.l l2002^ . We 
use bimodal distribution with equal fraction of objects in 
each mode. An average velocity in the first mode is 175 
km s ~^ and in the second it is 750 km s distribution 
in each mode is maxwellian. 

For the scenario without kick we proceed as follows. 
For the second group of tracks we used ranges indicated 
in the table 1. Width given in the table is calculated as 
0.5(max-min)/(max+min). For the first family of tracks 
we used the range for a from 230 to 346 Rq , for Mi from 
8.4 to 12.6 Mq, and for M2 from 1.6 to 2.4 Mq. 

For the scenario which takes into account an additional 
velocity gained by a NS at birth we used wider range of 
initial semimajor axis: from 200 to 2000 Rq. Masses are 
chosen in the same way as for the zero kick variant. 

The results of the calculations for non-zero kick are the 
following (we assume the total number of all NSs in the 
Galaxy as 10^, and the galactic age as 1.5 10^° yrs). In the 
first channel (fig. 1 left panel) we do not obtain significant 
number of massive NSs. Most of these objects are formed 
in the second channel. Formation rate of massive NSs was 
found to be 6.7 10~^ yrs"^. This corresponds to ^ 10 000 
of these compact stars in the Galaxy. For zero kick the 
formation rate is larger 4 10~^ yrs"^, so the total number 
is - 60 000. 

Certainly only a fraction of massive NSs at any given 
moment passes through stages which are observable, ie. 
the accretor stage and the stage of radio pulsar. Some 
of these objects are at stages of ejector and propeller or 
georotator. All three of them are not favourable for detec- 
tion^. In the table 2 we give fractions of massive NSs on 
each stage. It is clear that accretors are more numerous 
(but the number of massive NSs at the stage of superEd- 
dington accretion is negligible). 

® We note, that the ejector stage does not coinside with the 
radio pulsar stage, but includes it as a substage. So here we 
are speaking about non-detectability of eje ctors which are not 
active as rad io pulsars. See for example llLipunovlll992l) or 
iLipunov et aLil996i') for more details. 



0,2 - 



0.15 - 




M(NS)/Mo 

Fig. 2. Mass distribution of NSs. As we are interested only 
in the massive population we do not show the results for 
compact objects with M < 1.8 Mq. Upper mass limit cor- 
responds to SkyS with maximum rotation (Ouyed 2004). 
The dashed line represents results for the scenario with 
zero kick. The solid line - non-zero kick. Left peaks for 
both distributions correspond to NSs with a single episode 
of accretion. Right peaks are formed by NSs which also in- 
creased their masses via accretion from WDs. Distribution 
were normalized to unity, ie. an area below each line is 
equal to one. 

For the non-zero kick model about 25% of accreting 
massive NSs have normal stars as secondaries, the rest 
75% have WD companions. For zero kick nearly all mas- 
sive NSs accrete from WDs which fill their Roche lobes. 

Mass distributions for both scenarios are shown in the 
fig. 2. Note, that all small details in the figure are due to 
statistical noise (for example, the first peak on the rising 
part of the dashed curve, or the middle peak on the solid 
one) . The only important details are the two peaks at M 
2.3 Mq and M ^ 3.1 Mq, which correspond to tracks 
2a and 2b (see the right panel of Fig. 1 and table 1). 
As we found only two groups of tracks which lead to the 
formation of massive NS only results obtained for these 
groups are shown. All contributions from other types of 
tracks are not considered here and in fig. 3. 

Finally, in the last figure we represent luminosity dis- 
tributions. For the scenario with non-zero kick about 
1/2 of massive NSs have M > 2.5 Mq. Taking all to- 
gether we can conclude that in the Galaxy there are sev- 
eral thousand of accreting massive NSs with luminosities 
10^4 ^ L < 1036 erg s^^. 

4. Discussion and additional comments 

Here at first we notice some uncertainties of the scenario. 
Then we briefly discuss a possibility of massive NS for- 
mation in globular clusters, low-mass BHs, and types of 
sources which can host massive NSs. 
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Fig. 3. Luminosity distribution of accreting massive NSs. 
The left bin includes all sources with L < lO'^" erg s~^. 
The dashed line corresponds to the scenario with zero kick. 
The solid line - non-zero kick. In the ranges 10^° ^ L ^ 
1034 ej-g g-i and 10^^ ^ L <, 10^^ erg s~^ the number of 
systems is not equal to zero, however, it is very small. All 
distributions are normalized to unity. 



4.1. Correlations between initial parameters of neutron 
stars 

The scenario of binary evolution that we use has different 
types of uncertainties. Here we touch just one of them - 
possible correlations between parameters of the scenario. 

In our calculations we assumed that such initial pa- 
rameters of NSs as spin period, magnetic field, mass, and 
velocity are uncorrelated with each other. The reason for 
this assumption is trivial: there is no direct indication on 
such correlations. However, theorists suggested plethora 
of them. We just give a hst (probably incomplete) of pos- 
sibilities and corresponding references to original papers. 



— Spin ~ magnetic field jThompson fc Duncanlll993l). 

— Magn etic field - mass (jPopov et al.ll2002l : iHeeer et alJ 
120041) . 

— Internal structure - yelocitv llBombaci fc PoDovl2004l) . 

— Binarity - velocity l|Podsiadlowski et alJ l2004j) . 

— Core mass - velocity fScheck et al. 2004|). 

As here we mainly speak about masses and kick ve- 
locities, let us make a short comment on last two items. 
Our calculations may not be strongly influenced by such 
correlations. The reason is that the mass added during ac- 
cretion is much larger then a difference in initial masses. 
I.e., as we need to accrete nearly two solar masses to ob- 
tain the most massive SkyS at maximum rotation, we can 
safely forget about the range of initial masses. That's why 
in our calculations we even assume, that all NSs have the 
same initial mass. For the same reason we can neglect the 
second item in the list above. 

As here we deal with systems with high accretion rates, 
and the magnetic field is assumed to decay due to accre- 



tion, the first item also is not very important for us. Even 
if a NS was a magnetar in its early years, later on after 
field decay it can follow a normal evolutionary path. 

We do not discuss here a possibility of phase transition 
in a NS interiors. Formation of a quark phase due to mass 
increase can significantly influence the following history of 
a binary system, but we just neglect it here as in that case 
formation of massive NSs is impossible. 

4.2. Globular clusters 

All evolutionary tracks that we present above correspond 
to binary evolution in the Galaxy, and so they cannot 
be directly applied to globular clusters. However, we are 
mainly interested in systems which after a formation of 
a NS appear to be stiff, ie. orbital velocities in binaries 
are larger than a velocity dispersion in a cluster. It is also 
true during the following evolution of a system, and it can 
be violated only at the stage of Roche lobe overflow by a 
WD. We can conclude that a dynamical influence of the 
globular cluster stellar population should not destroy sys- 
tems under discussion. Still duration of various evolution- 
ary stages can be different in clusters and in the galactic 
disc, and so our estimates of relative fractions cannot be 
valid for globular clusters. 

It is possible to speculate that as the formation rate of 
millisecond pulsars is enhanced in globular clusters then 
the formation rate of massive NSs can also be higher there 
than in the disc of the Galaxy. It is an important question 
because massive NSs from globular clusters can enrich the 
disc population of these objects (so, our calculations which 
neglect this contribution give a lower limit on the number 
of massive NSs in the Galaxy). In our opinion millisecond 
radio pulsars and X-ray sources in globular clusters can 
be good candidates for a search of massive NSs. 

4.3. Low-mass black holes 

As it was described in the Introduction at present all well 
determined values of NS masses lie below ~1.5 Mq, on 
the other hand most of BH mass determinations lie around 
values 6-10 Mq (Ziolkowski 2004). So there is an indica- 
tion on a gap in intermediate mass range. Briefly we can 
say, that accretion cannot fill this gap if, as it is standardly 
assumed, NSs are formed with M ^ 2 Mq and BHs are 
formed with M <; 5 Mq . 

If an EOS of NSs with a very high Mmax is realized in 
nature, then up to ~ 3 Mq or even further, in the case of 
maximum rotation, we can find NSs. Otherwise the gap 
above ^ 2Mq should be populated only by BHs. Even in 
the case of the EOS discussed by Ou ved fc B utler ( 1 999[l 
low-mass BHs can form from rapidly rotating massive NSs 
as they slow down. 

Fig. 2 (the solid line) clearly shows that the number of 
low-mass BHs (or any other type of compact objects) with 
M ^ 3.2 in our scenarion is small. However, if A/max is ~ 
2 Mq, and if a binary is not significantly influenced during 
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a BH formation (i.e. accretion continues), then the number 
of BHs with Mmax ^ M ^ 3.2 Mq is non-neghgible."^ 

There are several examples of binary systems with 
an estimate of a mass of a cornpact object ~ 3-4 Mq 
I Orosz et al.l2nn4 : Shahb az et alJ20n4l) . These objects are 
considered as BH candidates. In principle such objects can 
be formed in the frame of the scenario discussed above af- 
ter a mass of a NS exceeds the Oppenheimer-Volkoff limit. 

4.4. Possible candidates 

Main astrophysical manifestations of massive NSs are the 
same as for normal NSs: X-ray sources and radio pulsars. 
However, there are differences. Very massive NSs should 
have short spin periods as they get an additional mass by 
accretion which spin-up them and provoke magnetic field 
decay. * Of course a given millisecond pulsar can contain 
a NS with a normal mass. Presence of a low-mass degen- 
erated companion (a WD) can be an indication that the 
system can hide a massive NS. One can mention another 
additional signature of a massive NS - very low magnetic 
field. 

If the magnetic field is very small, then the Alfven 
radius becomes less than the NS radius, and the accretion 
disk can nearly approach the NS surface. This situation 
takes place when 

here M_8 = M/10~^ Mg/yr, rio = tns/IO km and m 
— mass of NS in Solar units. Thus, the magnetic field 
strenght can be ^ 10^ G for Eddington accretion rate. In 
that case a formation of a boundary layer is favorable, and 
in the NS s pectrum an additional thermal component can 
be present l)lnogamov fc Sunvaevl [l999). For massive NSs 
including Ouyed EOS radius of a star is smaller than the 
distance to the last stable orbit, so the disc cannot actually 
smoothly approach the surface but qualitative properties 
of a spectrum will remain the same. 

All these consideration can be summarized in a list of 
types of objects which can contain massive NSs. 

— X-ray sources with weak pulsations with signatures of 
a boundary layer; 

— millisecond X-ray pulsars with WD companions; 

— millisecond radio pulsars with WD companions; 

— other kinds of millisecond X-ray pulsars; 

— other kinds of millisecond radio pulsars. 

By "other kinds" we mean millisecond pulsars with other 
types of companions or isolated (but old) ones. Note, 
that we do not include into our calculations secondary 

Mass growth of NSs and BHs in close binaries is also dis- 
cussed in iBoeomazov et aLllioOsI) and in (Bogomazov et al., 
in preparation). 

* If a NS has a very short spin period then pulsations in an 
X-ray source can be undetectable as it is observed for many 
low-mass X-ray binaries. 



companions with very low initial mass (brown dwarfs). 
However, such systems cannot produce massive rapidly 
rotating SkyS as the total amount of accreted matter is 
not sufficient. NSs with very low-mass companions like 
the millisecond accreting pulsar SAX J1808. 4-3658 or like 
"black widow" -like radio pulsars can be produced in our 
scenario via evaporating degenerated or non-degenerated 
se condaries (see discussions on the evolution of this source 
in lErgma fc AntiDOvall999tlBildsten fc Chakrabartvl20nil 
and references therein). 

Unfortunately our calculations cannot provide exact 
numbers of objects of each type. Uncertainties are con- 
nected with influence of population of sources from glob- 
ular clusters and with uncertainties of the scenario itself. 
For example we absolutely do not take into account in- 
fluence of rotation on the evolution of normal stars (see 
iLanger et ablEnni. 

lOuvedl (j2002L l2004^ discussed three binary systems as 
possible candidates to massive SkyS: 4U 0614H-09, 4U 
1636-53, 4U 1820-30. From the point of view of evolu- 
tionary scenarios discussed above all three really can con- 
tain a massive NS. 4U 1820-30 is especially interesting. 
The orbital period of the system is only 11 minutes which 
means that the secondary is a l ow- mass helium star (see 
iBallantvne fc Strohmaven l2004l and references therein) . 
However, this sources is situated in a globular cluster, and 
so our considerations should be applied with care. 

5. Conclusions 

We discussed possible channels of massive NS formation. 
If the EOS bas e d on the Skyrme model suggested by 
lOiived fc Butleil l)l999l) is reahzed in nature then these 
objects can be SkyS with masses up to 3.45 Mq for max- 
imum rotation. The estimated numbers of these sources 
in the Galaxy is high enough. Most favourable candidates 
are X-ray binaries with WDs as donors, millisecond radio 
pulsars in pair with WDs and accreting NSs with very low 
estimated magnetic field. If no of so massive NSs are found 
in these systems then the SkyS EOS has to be rejected. 
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